Two-photon Compton process in pulsed intense laser fields 
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Based on strong-field QED in the Furry picture we use the Dirac-Volkov propagator to derive 
a compact expression for the differential emission probability of the two-photon Compton process 
in a pulsed intense laser field. The relation of real and virtual intermediate states is discussed, 
and the natural regularization of the on-shell contributions due to the finite laser pulse is high- 
lighted. The inclusive two-photon spectrum is two orders of magnitude stronger than expected from 
a perturbative estimate. 
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Recently, experiments on two-photon emission by elec- 
trons in an intense laser field have been proposed [T] [2] . 
The motivation of such experiments is seen in [TJ [3J in an 
attempt to verify the Unruh radiation [5] which is re- 
lated to the physical vacuum experienced by accelerated 
observers in a flat space-time and manifests as the emis- 
sion of entangled photon pairs off accelerated charges. 
The relation between the QED two-photon Compton pro- 
cess and Unruh radiation has been analyzed further in 
[6]. As the intensity of the Unruh radiation increases 
with the acceleration [7], e.g. caused by a strong electro- 
magnetic field, one is naturally lead to think on strong 
laser fields as driving force. High field strengths are nowa- 
days achieved by the temporal compression of laser pulses 
to a few femtoseconds. The temporal structure of the 
laser pulse has a strong impact on the nonlinear Comp- 
ton spectrum [5] (hereafter termed one-photon Compton 
scattering) and its cross channels [§] [ID] • 

The theoretical description of the two-photon Comp- 
ton process (also termed double Compton scattering) 
in the perturbative (weak-field) regime had been ac- 
complished first in 112) and was verified experimen- 
tally soon afterwards |13j : for a more recent experiment 
cf. [H]. 

In second-order strong-field processes (e.g. [T31tT5] h in- 
termediate particles can become real (i.e. go on their 
mass shell) due to the presence of the background field. 
The on-shell contributions have been discussed [HI [20] as 
Oleinik resonance singularities. The relation of the on- 
and off-shell processes in the case of a photon propaga- 
tor has been analyzed recently in [T5J [21] for the trident 
pair production with respect to the interpretation of the 
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FIG. 1: Feynman diagrams for the two-photon Compton pro- 
cess, where two photons (wavy lines) are emitted. The double 
lines representing laser dressed Volkov in and out states and 
the Dirac Volkov propagator between the two vertices, respec- 
tively. 



SLAC E-144 data [22]. 

In this paper, we provide a complete description of 
the nonperturbative two-photon Compton process in a 
pulsed intense laser field. We show the significant mod- 
ification of the two-photon Compton process by short 
laser pulses when compared to infinite plane-wave fields 
considered previously [T5]. For the first time, an exact 
description of this process involving a temporally shaped 
laser-dressed Dirac-Volkov propagator within the Furry 
picture is given. We calculate the photon pair emission 
probability and compare it to the one-photon Compton 
probability, finding a substantially increased two-photon 
yield in short intense laser pulses as compared to pertur- 
bative estimates. 

The Feynman diagrams in Fig. [T] correspond to the S 
matrix 

S = -ze 2 J d 4 xd 4 y%, {y)f 2 ^Lg{y, x)^^l%(x) 
+ (1**2), (1) 

where (1 ■<-» 2) means exchange of photons 1 and 2 ac- 
counting for the symmetrization of the two-photon wave 
function. ty p ("ly) is the Volkov wave function for an 
electron in the entrance (exit) channel with momentum 
p (p'), and Q(y,x) denotes the dressed electron prop- 
agator [23], ki (ej) stand for the four-momenta (four- 
polarizations) of the emitted photons. A dot, e.g. in ki-x, 
denotes the scalar product of four- vectors, and ^. = 7-ej. 

We describe a linearly polarized laser pulse by the four- 
potential A p — E e^a((j))/(jj with polarization four- vector 
f} 1 = (0,1,0,0) and four-momentum fc M = (w, 0, 0, — ui) 
with e • k = 0, phase cj) = k • x, a(4>) = g(4>) cos</> and 
temporal envelope function g((f>) which vanishes as \(f>\ — ¥ 
oo. The nonlinearity parameter is clq = eEo/mui with 
laser frequency u, electron mass to, charge e, and peak 
electric field Eq. The perturbative regime corresponds to 
a < 1. 

Both, ^ p (x) and Q(x,y), depend on the Volkov ma- 
trix functions E p (x) = r p (</>) exp{— ip ■ x — if p {4>)} with 
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r p = 1 + e$4/ (2k ■ p) and the nonlinear phase 



f p (<b) = a p / #'a(^)+^ / #'aV) (2) 



with coefficients a p — ma (p ■ e)/(p ■ k) and /3 p = 
m 2 a^/(2p • k). We employ light cone coordinates with 
p± = p° ±p 3 , p^ = (p 1 ,p 2 ) and p = (p+,P_l) such that 
fc_ is the only non-vanishing component of the laser four- 
momentum with (f> — k-x + /2 and d 4 x — k^dtfidx-dx.^. 

The integrations over the transverse and — components 
provide momentum conserving delta distributions, such 
that the S matrix ([!]) becomes 
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xT Pl ('<p){* 1 r p (ip)u p + (It* 2), (3) 

with Pi = p — k\ and s = (p'_ + k\- + k,2- — P- ) /fc- > 0. 
For certain pulse profiles a part of the non-linear phases 
f p referring to the ponderomotive energy might be rewrit- 
ten to generate explicitly a mass contribution in the de- 
nominator, leading to the familiar form of the propagator 
with a mass-shift Am 2 — m 2 a§/2 in infinite plane waves, 
e.g. in [15]. For our purposes, ^ is more suitable. The 
variables s and I are continuous analogues to the number 
of exchanged photons. The integral over I in ^ accu- 
mulates all possible paths the system may take. It can 
be evaluated analytically by the contour integration tech- 
nique: The integrand has a single pole at £ — l\ — ie' with 
l x = (m 2 - P 2 ) /2k ■ P x and e' = esign(P 1+ ). To safely 
apply the residue theorem, one has to transform the inte- 
grand according to £/(£ - 4 + is') 1 +£i/(£ - £\ + is') 
such that the nontrivial part goes to zero as £ — > oo. The 
result of this procedure, 
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takes the structure of a fermion propagator in light-front 
form [21] due to the integration over the light-cone com- 
ponent fc_. In the second line of Q a time ordering 
(in the laser phase) is introduced in the sense that the 
emission at the second vertex has to happen at a later 
time than the emission at the first vertex by means of 
the step function 8((f> — ifj). The quantity £\ controls 
the amount of momentum £\k transferred to the electron 
such that the propagator momentum is on its mass shell, 
i.e. (Pi + £\k) 2 = m 2 . Additionally Q includes a part 
oc 6(4> — if>) where both photons are emitted simultane- 
ously. This " instantaneous propagator" , which is also 



known as light-front zero-mode propagator, is specific to 
the fermion propagator and does not appear in the analy- 
sis of the trident process [21] with the photon propagator. 
The propagator pole in Q always lies below the real axis 
due to momentum conservation. In particular, we have 
p' + = P\ + — k 2 + > with fc 2+ > and therefore P i+ > 0. 
Negative values of P + would shift the pole to the upper 
half plane corresponding to the opposite time ordering 

The six-fold differential probability of two-photon 
emission per incident laser pulse finally reads 
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with a = e 2 /4ir and the phase integrals 

A n (s) = J d(/>a n ((l>)e i3 ' t ^ if * w+if >'' w , (7) 
B nl (s,£ 1 )= J ##0(<£-V>K(£)a i (V') (8) 

and the Dirac structures To = /* 2 ^./\, Pi = Xpi/* 2 ^/* x + 
i\U\X p , P 2 _= 4d p d p ,(e* 2 ■ k){e\ ■ k% Poo = £Gi£, 
^oi = f 2 G\(X p 1 /* 1 + /* 1 X p ), Uio = (Xptfe+tlXpJGi^l, 

Ull = ( X p'f* 2 + f2 X Pi)Gl{X Pl f 1 + f^Xp), U 02 = 

f 2 G- i _Xp 1 f 1 X pi U20 — X p >/* 2 X Pl Gi/* 1 , U12 = [X p 'f* 2 + 
^Xp 1 )GiXp 1 ^*X p , U21 = X p i f 2 X Pl G\(X Pl f 1 -\- f[Xp) , 
U22 — 8dpd p >dp i (e 2 -k)(el-k)(Pi-k)k', where the abbrevia- 
tions X p = d p l/:/, dp = aom/ (2k-p) and Gi = rf > 1 +£ik'+m 
are employed. The phase integrals Aq, Bqi, B n o and Poo 
are numerically non-convergent because of the missing 
pre-exponential factor. However, these integrals can 
be defined as a superposition of convergent phase in- 
tegrals by applying a quantum gauge transformation 
ei -> ei + Xih [5T], yielding, e.g. (s - £i)B ol (s,£x) = 
iAi{s) + (a Pl - a p /)Pii(s,4) + {fin - M B 2i(s,£i)- 
Thus, gauge invariance reduces the number of indepen- 
dent phase integrals from twelve to six well-behaved ones. 

Applying the Sokhotsky-Weierstrass theorem, (I— £\ + 
ie)- 1 =V(£-£ 1 )- 1 ~i7r5(£-£ 1 ), in @ one identifies the 
imaginary part of ^# as the contribution of real (on-shell) 
electrons with the aid of the delta distribution. By recall- 
ing that principal value integration (P) effectively means 
cutting out a small interval around zero it becomes clear 
that the real part of M refers to the virtual (off-shell) 
process. The phase integrals A n contribute only to the 
off-shell process. Employing the completeness relation 



for spinors G\ = Up^e^up^i^a to the numera- 
tor of the propagator (a is the spin of the intermediate 
electron), the on-shell part of S yields a factorization on 
the amplitude level 
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where denotes the S matrix for the one-photon 

Compton process with the given energy-momentum con- 
servation (cf. equation (29) in [8]) integrated over all in- 
termediate states. Numerically, the interference between 
the two Feynman diagrams in Fig.[T]turns out to be of the 
same order of magnitude as the non-interference terms. 

In the weak-field limit ao <C 1, the off-shell part of the 
S matrix is proportional to ao in lowest order correspond- 
ing to the absorption of one laser photon reproducing the 
known perturbative result. Contrarily, the leading order 
contribution to the on-shell part is oc Oq, since at least 
two laser photons are needed to satisfy the on-shell en- 
ergy momentum conservation. In the limit of infinite 
laser waves the on-shell part of the two-photon rate gives 
rise to Oleinik resonances [19] at frequencies 
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where rii — (1, cos ipi sin 9i, simpi sin 9i, cos 9i) is the unit 
vector in the direction of ki (with 9i and ipi denoting the 
usual polar and azimuthal angles of the emitted photons 
and integer £). The resonances for 07 (W2) emerge when 
the propagator of the first (second) Feynman diagram in 
Fig. [T] comes on its mass shell. A detailed discussion of 
the resonance behavior will be given elsewhere; it has its 
own right but does not affect the following results. 

For our numerical evaluations we consider electrons 
with a Lorentz factor 7 = p Q /m = 10 4 , available, e.g. at 
the European XFEL electron beam [25], in head-on col- 
lisions with the laser pulse. Calculations have been per- 
formed for ao = 1 and a pulse shape g(<fi) — cos 2 (^) 
for — t > cf> > t and zero otherwise, such that r is the 
dimensionless FWHM pulse length with r = 20 corre- 
sponding to 9 fs FWHM for ui — 1.55 eV. (A similar 
kinematic situation with the same center-of-mass energy 
could be achieved by colliding an XFEL X-ray pulse [25] 
with low energy electrons, e.g. 7 = 10, provided by an op- 
tical laser acceleration set-up [26].) In Fig.[2j results are 
exhibited for the differential probability for two-photon 
emission as a function of wi and 101. Since the motion 
of the electron is relativistic, the radiation is produced 
in a cone around the spatial direction of p with a typical 
opening angle of I/7. We show the differential probabil- 
ity at angles #1,2 = 1/7 and ipi = tt/2 and </?2 = Stt/2, 
i.e. the two photons are emitted in a plane perpendicular 
to the polarization plane of the laser. The left panel in 




0.4 0.6 0.8 1.0 0.2 0.4 0.6 

Ui [GeV] uii [GeV] 

FIG. 2: Intensity distribution for two-photon emission in the 
ui\-u)2 plane. We show the complete emission probability [a)] 
and the off-shell contribution [b)]. The color code represents 
the logarithm of the six-fold differential probability in eV, 
log 10 (d 6 W / ' du>idfliduj2dfl2 [eV]). For parameters see the text. 

Fig. [2] shows the complete differential probability as the 
sum of on- and off-shell parts, while in the right panel 
the off-shell part is exhibited. The probability distribu- 
tions display complex characteristic patterns. The dif- 
ferential spectrum in Fig. [2] a) is dominated by the on- 
shell part in almost the whole W1-W2 phase space for both 
uji > 200 MeV, where it is roughly one order of mag- 
nitude larger than the off-shell part. This is a generic 
feature also for different scattering angles. The on-shell 
part shows a rectangular pattern, which is aligned par- 
allel to the coordinate axes. The spectrum has maxima 



in regions where the Oleinik resonances ( 10 ) would occur 



for infinite plane-waves and is particularly strong where 
both types of resonances intersect. The resonances (10) 
for £ = 1 are indicated in Fig. [2] a) as dotted lines. The 
off-shell part exceeds the on-shell part for at least one of 
the u>i below 200 MeV, where the maxima of the distri- 
bution are roughly aligned with the different harmonics 
£ of the infinite plane- wave energy correlation which read 



u-p-p-ki-(i+ig) 



'k ■ ki 



Ak-p 



k-h 



(11) 
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these are shown as dotted lines for I = 1, 2 in Fig. [2]b) . 
For higher photon energies the pattern is more irregular 
but still symmetric with respect to an exchange of u>i and 
u! 2 - Fig.[2]evidences the striking differences to the infinite 
plane- wave case: The strict W2(wi) correlation of (11) 
gets completely lost. Instead, resonance like structures 
with subpeaks appear which are produced by the pon- 
deromotive broadening mechanism, resembling the ones 
observed in the one-photon Compton process [8]. 

In Fig. [3] we exhibit the inclusive spectrum 
(PW/duJidfli arising from ^ by integrating over the 
phase space of photon 2. For soft photons 0J2 — > the 
emission probability of two-photon emission becomes di- 
vergent. The cancellation of this infrared divergence by 
soft virtual photons due to loop corrections of one-photon 
scattering is ensured by the Bloch-Nordsieck theorem [57] 
as in the perturbative case. For practical purposes, how- 
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ever, we include an infrared cutoff w™ m = 100 keV to 
avoid the soft-photon divergence in the spirit of [15] . The 
value of the integral is rather insensitive to a variation of 
the cutoff in the range of 1 — 1000 keV. 

This inclusive spectrum accounts for the experimen- 
tal observation of only one of the two photons. To 
compare with one-photon Compton backscattering, we 
choose 6\ = (p\ = 0. In the case of strong laser fields, 
e.g. for arj = 1, the inclusive spectrum is found about two 
orders of magnitude below the one-photon spectrum for 
lji > 200 MeV. At photon energies wi < 200 MeV, the 
two-photon process exceeds the one-photon process (see 
Fig. [3|, opening, at least in principle, a window to access 
its observation without coincidence measurements. 

Approximating = fd^^ » £ ^1 

and integrating over u>\ , one can estimate the total num- 
ber of produced pairs as 1.1 x 10 -3 per pulse and electron 
as compared to 5 x 10~ 2 coming from the single photon 
process. With an assumed laser repetition rate of 10 Hz 
one can expect 950 two-photon events as compared to 
43000 single photon events in one day which should be 
sufficient for an experimental observation. The coinci- 
dence detection of rare two-photon events, where both 
photons are emitted within a small opening angle has 
been successfully demonstrated in the photon splitting 
process [28] ten years ago. The experimental sensitivity 
might be increased by a simultaneous detection of the 
scattered electrons (like in photon tagging). The electron 
beam should be a dilute beam tuned to one interaction 
per laser pulse. 

In the weak-field regime the rate of the two-photon 
Compton process is suppressed by a factor of a(k-p/m 2 ) 2 
relative to the one-photon Compton process for k-p <C m 2 
1 1 1 j . as for our kinematics. For the momenta considered 
here the estimated suppression in the weak field regime 
is 3 x 10 -5 . To discuss the relevance of the two-photon 
emission in strong laser fields we define the two- to one- 
photon ratio as TZ = For r = 10 we 
obtain a value of TZ = 10~ 2 at ao = 1 which is about 
two orders of magnitude larger than the perturbative es- 
timate. For lower values of a < 0.1 the suppression of 
the two-photon probability rapidly approaches a constant 
value of TZ = 10~ 4 as anticipated in [11]. Considering the 
different contributions we find that the ratio for the on- 
shell process TZ on = 0.01a 2 , for ao < 1; the on-shell ratio 
is independent of ao for ao < 0.1, and above ao > 0.1 the 
value increases and reaches TZ g — 10 -3 at ao = 1. 

Our approach furthermore opens the avenue towards 
a detailed study of the two-photon polarization which is 
considered as a signature of the Unruh effect in [TJ [3]. 
The two-photon emission as a QED process in itself is 
also interesting with respect to the quantum radiation 
reaction [25] as multiple incoherent one-photon Compton 
scatterings. The sequential Compton scattering appears 
thereby as a factorization of the resonant on-shell part of 
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FIG. 3: Inclusive spectrum vs. the one-photon Compton spec- 
trum as a function of the emitted photon frequency wi. 

the scattering matrix element on the amplitude level and 
is complemented by the possibility of coherent emission 
due to the off-shell part. Our numerical calculations show 
that up to 30 % of the total photons are due to the off- 
shell process and therefore beyond a description based 
solely on real intermediate photons. 

In summary we provide the first complete evaluation 
of the differential probability of two-photon emission by 
an electron in a short intense laser pulse. The on-shell 
part of the matrix element factorizes into subsequent 
one-photon Compton processes and gives naturally a fi- 
nite contribution to the differential probability due to the 
temporal pulse structure. Our result allows for the first 
time an unambiguous comparison of the probability of 
the two-photon process in relation to the one-photon pro- 
cess. We find an increased two-photon yield by two orders 
of magnitude as compared to the perturbative estimate 
of two-photon Compton scattering even for moderately 
strong laser fields ao ~ 1 which are available presently at 
various laser facilities. 
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